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I present quenched domain wall fermion and 2+1 flavor improved Kogut-Susskind fermion calculations of the 
hadronic vacuum polarization which are used to calculate the 0(o?) hadronic contribution to the anomalous 
magnetic moment of the muon. Together with previous quenched calcuations, the new results confirm that in the 
quenched theory the hadronic contribution is signifcantly smaller (~ 30%) than the value obtained from the total 
cross section of e + e~ annhilation to hadrons. The 2+1 flavor results show an increasing contribution to g — 2 as 
the quark mass is reduced. 



1. Introduction 

The anomalous magnetic moment of the muon 
(g-2) has been measured to incredible accuracy at 
Brookhaven National Lab's E861 experiment 1 
and has also been calculated precisely in the Stan- 
dard Model using the dispersion relation for the 
hadronic vacuum polarization, n(</ 2 ), and the to- 
tal cross section for e + e~ annihilation to hadrons 
(isospin symmetry allows the decay rate for r lep- 
tons to be used as well) 2 3 . The theory calcu- 
lations agree with experiment at roughly 2.7 and 
1.4 standard deviations, depending on whether 
solely e + e~ data are used or if r decay data are 
included as well. 

A while ago it was shown how to compute 
the lowest order (in a) hadronic contribution 
to g — 2 completely from first principles, using 
lattice gauge theory to calculate the hadronic 
vacuum polarization in conjunction with contin- 
uum perturbation theory in the QED coupling 
constant 13] ■ In this work we continue along this 
path, past new high statistics quenched domain 
wall fermion calculations, to 2+1 flavor improved 
Kogut-Susskind fermion calculations on large vol- 
umes V 5s (3.5 fm) 3 and with small quark masses 
[mi = (m u +TTid) /2 ss m s /10). In the latter case I 
use configurations from the MILC collaboration. 
The 2+1 flavor calculation is not complete, so I 
do not give numbers for g-2 in this case. After 
presenting results, I discuss the work that still 
needs doing and prospects for a precision lattice 
calculation of g-2. 

The vacuum polarization is defined from the 
two-point correlation function of the electro- 



magnetic current 
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where q is photon momentum. See |4I6| for de- 
tails of the lattice calculation. The central idea 
is to calculate n(g 2 ) on the lattice, and fit it to 
obtain a continuous function that can be read- 
ily integrated (numerically) with a know function 
from continuum perturbation theory from q 2 = 
up to a cut-off ~ (1/a) 2 , the rest of the integral 
being done using perturbation theory. For the 
quenched case, the (continuum) form of n(g 2 ) is 
known and provides an ansatz for the fit[Hj. 

n(q 2 ) = —ft+C\n(a 2 (q 2 +f, 2 )), (2) 
q* + my 

where the first term comes from the vector meson 
bound state (delta function) and the second from 
a two particle continuum (cut, q 2 > /i 2 ). Note 
that fy and my can be extracted in the usual 
way from the zero momentum correlator as was 
done in |S] which leads to a more accurate fit, or 
treated as free parameters (here, I treat fv as a 
free parameter and take my from[Jj). This ansatz 
fits the data well and leads to (statistically) accu- 
rate results that are extrapolated to q 2 = 0, which 
is important because the low q 2 region dominates 
the one-loop integral that gives the hadronic con- 
tribution. For the 2+1 flavor case there is, of 
course, no such phyical ansatz, excluding the ex- 
perimental one which is, after all, the thing we're 
trying to compute in the first place. In this case 
I try several forms: the pole fit just described, 
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a simple polynomial, a log for the low q 2 region, 
and combinations of these. 

2. Results and Discussion 

The quenched calculations with domain wall 
fermions were done with inverse spacing a" 1 ~ 
1.3 and 2 GeV on 16 3 x 32 lattices, and m va i = 
0.02 and 0.04. H(q 2 ) is shown in Figure Q and 
values of = (g — 2)/2 for a single quark with 
unit charge are sumarized in Table 1. After in- 
cluding the u, d, and s quark charges, the total 
contribution for three degenerate quarks is con- 
sistent with previous quenched results |4"l6*j . con- 
firming that the hadronic contribution to g — 2 
is significantly less (~ 30%) than in the real 
world[2]. Thus it is probably not worthwhile to 
pursue further quenched calculations designed to 
eliminate systematic errors due to finite volume, 
non-zero lattice spacing, and unphysically large 
quark mass (which has been done to some extent 
in |4I6|L 
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Figure 1. The hadronic vacuum polarization for 
quenched domain wall fermions, a^ 1 = 1.3 GeV, 
mvai = 0.02 (circles) and 0.04 (squares). The line 
shows the result from continuum 3-loop perturba- 
tion theoryjH] for comparison. 



Instead, I have started a calculation on the 
2+1 flavor lattices from the MILC collaboration 



that were generated using a 2 — tad fermions (here 
I omit the Naik term for the valence quarks 5 ). 
Initial results were presented last yearj^l, and cur- 
rently II(g 2 ) is being calculated on new 40 3 x 96 
lattices with m; = 0.0031. Of course, these are 
very aggressive parameters, so the lattice genera- 
tion is somewhat slow. At the time of the meet- 
ing some 84 configurations existed on which n(<7 2 ) 
was calculated from a point-split current from a 
single site (and its nearest neighbors). These were 
separated by six trajectories and so are probably 
not independent. To improve statistics, I have 
begun calculating on a time-slice that is one-half 
the lattice size distant from the first. In addition 
the MILC collaboration plans to at least triple 
the length of the evolution (~ 3000 trajectories). 

In Figure 2 II(g 2 ) is shown for m va i = m; = 
0.0031 (40 3 x 96) and 0.0062 (28 3 x 96), or 0.1 and 
0.2 times the strange quark mass, respectively. 
Note that the large volumes and time sizes used 
here lead to very small values of q 2 which is quite 
important. The data show a slight increase as 
q 2 — > as mi decreases by a factor of two. While 
it appears small, a slight increase in this region 
changes the contribution to g — 2 significantly, so 
the fit in this region must be quite accurate and 
precise. In Figure 3 I show covariant polynomial 
fits to n(g 2 ) which tend to under-predict the data 
as q 2 — > (the x 2 of the cubic and quartic fits is 
acceptable, though). Pole fits like Eq. 2 do a bit 
worse, and a log fit does about the same. It seems 
that a better fit ansatz is needed. 

Chiral perturbation theory offers the means to 
understand the mass dependence and the small q 2 
behavior of H(q 2 )- This result does not appear in 
the literature, to the best of my knowledge, and 
even if it did, to be effective in this case, it would 
probably have to be augmented with so-called 
staggered chiral perturbation theory (SxPT) [§] 
to obtain accurate results. Indeed, recently the 
MILC collaboration has performed a fit to all of 
their 2+1 flavor data for the pseudo-scalar decay 
constant using SxPT^Uj. The combination of a 
large data set and S%PT allowed for a ~ 1-2% de- 
termination of Jk/ f-n which is quite remarkable. 
Since the vacuum polarization can be measured 
accurately on these same lattices, perhaps a sim- 
ilar fit will prove to be as precise here as there. I 
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am now considering such an analysis for II (q 2 ). 



Table 1 

The magnetic anomaly et M for a single flavor with 
unit charge from quenched domain wall fermions. 



a- 1 


m va i 


configs 


(a?) 




1.31 


0.02 


337 


750(35) xlO" 


-10 


1.31 


0.04 


337 


669(23) xlO" 


-10 


1.98 


0.02 


299 


730(51) xlQ" 


-10 
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Figure 2. Same as Figure 1 except for 2+1 flavor 
improved Kogut-Susskind fermions, a -1 rj 2.3 
GeV, m s = 0.031, m va i = mi = 0.0031 (open 
circles) and 0.0062 (filled squares). 
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Figure 3. Polynomial fits to U(q 2 ) for the 2+1 
flavor case with m va i = 0.0062. 
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